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Abstract 



We analyze the perturbative expansion of a condensate in the 0{N) non- linear sigma 
model for large N on a two dimensional finite lattice. On an infinite volume this 
expansion is affected by an infrared renormalon. We extrapolate this analysis to the 
case of the gluon condensate of Yang- Mills theory and argue that infrared renormalons 
can be detected by performing perturbative studies even on relatively small lattices. 
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1 Introduction 



In SU(3) Yang-Mills theory on the lattice, the coefficients of the perturbative expansion of the 
simplest observables, the elementary Wilson plaquettes, have been computed to large order ]l], [2|] 
and they show evidence of factorial growth. These perturbative calculations are needed not only 
to clarify the presence and relevance of infrared (IR) renormalons Q but also to extract the values 
of physical quantities, such as heavy quark decay constants Q from the results of Monte Carlo 
simulations. 

The perturbative coefficients of the elementary Wilson plaquette, an observable corresponding 
to the SU(3) gluon condensate, are known analytically J|] up to the third order in as, the weak- 
coupling parameter at the ultraviolet lattice scale, related to the usual lattice constant (3 by Anas = 
6//?. At larger orders these coefficients have been obtained by a numerical method based on the 
stochastic formulation of the theory a la Parisi-Wu 0]. One considers the weak coupling expansion 
of the Langevin equations and solves numerically the truncated set of equations corresponding to a 
given perturbative order. The coefficients have been computed numerically up to the order a s on a 
8 lattice and up to the order a s on a 12 4 lattice. In the sequel we shall always consider symmetric 
square (d=2) or hypercubic (d=4) lattices with M lattice points in each direction. 

Up to the sixth order the differences between the coefficients for M = 8 and M = 12 are found 
to be contained within the statistical errors, at present smaller than 10% Q 

An important question in lattice calculation is to what extent the finite volume affects the 
results. In particular, if the perturbative part of an observable is affected by an IR renormalon, 
the factorial behaviour of its coefficients is tamed by the finite volume. Then we should face the 
question| whether the lattice size (M = 8, 12) considered in p are sufficiently large to detect the 
IR renormalon factorial growth. 

In order to have a simple estimation of the finite volume effect we consider the typical integral 
involving an IR renormalon and we introduce an infrared cutoff related to the finite volume. For a 
condensate with dimension 5 we consider 

f klv dk 2 ( k \ 5 /l2 , , vr , 2tt 

with k nv and k ir the infrared and ultraviolet cutoff on the lattice of size M and spacing a. Using 
the one loop form of the running coupling 

as{k2) = l + b asHkykl v ) ' 

one obtains the formal expansion 

Ws(* S ,M) = £4 c e (M) , c e (M) = J* ^ (A)' tf-i 0^ . ( 2 ) 

The factorial growth of the coefficients is obtained from values of k around the steepest descent 
point 

^sd = ^uv e ^ • (3) 

t Since the stochastic method requires the introduction of two copies of the field for each order in as, the amount 
of computer memory grows proportionally to M 4 x number of loops. This is why only rather small lattices have been 
considered so far. 

* This problem has been stressed to us by Martin Beneke. 
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The IR renormalon factorial growth is then obtained only if fc S( j > k{ T , i.e. the steepest descent 
point lies inside the integration region, which gives the upper bound 

1 < 4d = $ In fy) + 1 . (4) 

Above this value the factorial growth is tamed. For M = 8, 12 and 5 = 4 one obtains £ s d — 6.5, 8.2. 
Therefore, among the coefficients computed in only the ones with i < 6 would be inside the 
bound. For higher order coefficients one should use lattices larger than M = 8. 

In this paper we try to probe the reliability of the indication (Q) and we shall argue that the 
factorial growth of the coefficients is still present at values of £ much larger than £ s d- 

First we consider the exactly solvable asymptotically free case given by the large N limit of 
the O(N) non- linear sigma model in a two dimensional finite lattice. We consider the perturbative 
expansion of the condensate with 5 = 2 and we compare the exact coefficients and their behaviour 
in M with the ones given by (|2|). We show that the large order coefficients agree with the ones 
obtained from the typical integral (Q) provided that the running coupling as(k 2 ) in the integrand is 
replaced by as(rk 2 ) with r < 1. We then apply this idea to the four dimensional case and consider 
the expression in (|l]) with 5 = 4. As suggested by the non-linear sigma model, we assume that 
the running coupling is taken at a scale smaller than k 2 . We then argue that the factorial growth 
extends to larger values of £ well inside the values considered in . 



2 Non- linear O(N) sigma model at large N 

IR renormalons are present in the perturbative expansion of any observable in the non-linear sigma 
model in two dimensions f6j. We consider the lattice version of the model | ona finite volume. 
The internal energy density U, which is a condensate of dimension 5 = 2, has the following operator 
product expansion 

U = a- 2 U pcvt + m 2 U! + m 4 U 2 + ■ ■ ■ , (5) 

with a the lattice UV cutoff and m the scale parameter generated by the breaking of scale invariance. 
In the limit N — > oo one has 

(am) 2 = 32 e" 4 ^ , (6) 

where (3 is the strong-coupling parameter at the lattice UV scale. The perturbative component has 
the following large iV expansion 

t/pert = 4vra + ^ W(a) + 0(N- 2 ) , (7) 

where we have defined the weak-coupling by Aira = 1/(3. 

The first 1/N correction contains IR renormalons. For finite volume, with M lattice size, and 
for infinite volume one finds 

r-x/a a 2 ( p k 47rQ , 



W(a, oo) 



it J a (2tt) 2 l + 4TTaA oc (k 1 ,k 2 ) ' 
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with 

. . , , . 1 x - k 2 - p 2 - (p + k) 2 



mi,m 2 r vx- . ■-/ ^ ^ 

1 ^ 2 ./-*/«, (2vr) 2 p 2 {p + k) 2 

As usual we have introduced the lattice momenta 

2 . a — , 2 . a. , , 

Pt = -sin-pi, p + fcj = - sin -(pi + fci), 

Ob Zi ( I _ 

with 

27T7Z ' 27T77Z ■ 

Pi = lff' h = l^d' n i ,m i = 0,±l,±2,---±(M/2-l),M/2. 

for M even and |pj| < 7r/a for M — ► oo. Notice that for finite volume the vanishing momenta do 
not contribute to the sums in (JsJ) and (||) and then one has an IR cutoff 

Pi, h < fe ir - — . 
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Figure 1: The ratio between A 2 §Q(ki, k 2 ) and its leading term of Eq. (11) plotted against ka; the 
spread for large k 2 is a lattice artifact representing the breaking of rotational invariance while the 
0{k 2 ) term accounts for the curvature. 



We now discuss the perturbative coefficients 



4tt 



W(a, Af) = $>' Q(M) , C e (M) = (-AnA M (k u k 2 )) t ~ 1 ■ 

Consider first the infinite volume case. For small momenta a\ki\ <C 1 one finds 

r 2 



Aoo(ki, k 2 ) 



7T" 

32 



(10) 



(11) 
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where k 2 = kf + k\ (see fig. 1). 

Therefore the integrand of (||) at small ki is 

4ira Ana 



1 + AnaAooik!, k 2 ) 1 + b a ln(rk 2 /k 2 v ) 



Ana(rk 2 ) , (12) 



with bo = 1 in this case (see ©). Therefore the integrand is given by the typical form (S) except 
that the running coupling is taken at the scale rk 2 smaller than k 2 . The presence of the Landau 
pole is directly responsible for the factorial behaviour at large order 

d(oo) ~ b e -\£- 1)1 r~\ £>1. (13) 



Non-leading corrections in l/£ can be estimated from (12) by taking into account that in the soft 
region the integrand A OQ (kik2) approaches the running coupling at a scale rk 2 rather than k 2 . 

We consider now the exact evaluation of the perturbative coefficients Ci(M) given by (0) both 
for M finite and infinite. 




Figure 2: The expansion coefficients Cg. for various lattice sizes M. 



In fig. 2 we report the exact coefficients Cg{M) as function of t for various M up to M = 256; 
for relatively large loops (£^4) and lattice size (M^16) the exact coefficients C((M) are fitted 
quite well by the expressions 




with 5 = 2 and r given by the soft limit ([Tl]) of the integrand r = 7r 2 /32 (see fig. 3). 

At finite M the factorial behaviour of the coefficients is tamed by the IR cutoff h ir . To see this 
we plot in fig. 4 the ratio of the exact coefficients Rn{M) = C 'e(M) / C t(oo) as a function of M for 
various I. As expected (see the discussion based on steepest descent estimation) this ratio vanishes 
by increasing M at fixed £ or by increasing I at fixed M. 
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Figure 3: The exact coefficients Ce and its approximate value given by Eq. (14) for M = 64 (+) 
and M = 256 (o). 



We come now to analyze the question whether a reliable estimate of the IR renormalon behaviour 
can be obtained from finite lattice calculations. In particular we analyze the limitation (||) suggested 
by the steepest descent estimation. 

In fig. 5 we plot £ s d as a function of M. At any given M we report also the maximum value of 
£ such that the exact coefficient Ce(M) is larger than a given fraction of the infinite volume exact 
coefficient C^(oo). We plot the maximum £ at which Ri(M) > a with a = 0.1, 0.2, and 0.5. 

The finite volume clearly tames the growth of the exact coefficients at large £. However the 
factorial growth does not stop abruptly and is still present for £ larger than the limiting value £ s ^ 
given in (||). For instance, for the values of M considered in fig. 5, one finds that the exact finite 
volume coefficients Ct(M) are larger than 20% of the infinite volume coefficient Cg(oo) for all values 
of £ with £&£ sd + 2. 



3 Extension to the Yang-Mills case and conclusions 

We consider a naive extension of the previous analysis to the gluon condensate in the Yang-Mills 



theory on a four dimensional lattice. To this end, taking into account that the expression (14) 
fits very well the exact coefficients in the non-linear sigma model, we assume that also the gluon 
condensate coefficients are approximated by the expression (0) with 5 = 4. 

For r we assume the same value as in the non-linear sigma model. In fig. 6 we plot the case of 
the gluon condensate. In the four dimensional expected, the effect of the infrared cutoff 

is less severe and the factorial growth is still present for much larger values of £. For instance, in 
the four dimensional case the finite volume coefficients are larger than 20% of the infinite volume 
coefficients for all values of £ with £ ^ £ sc j + 4. 
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Figure 4: The ratios Q(M)/Q(oo). 



We conclude that, if our fitted expression (14) for the gluon condensate perturbative coefficients 



is correct, the eight loop coefficient numerically computed in Q is of the order of 50% of the infinite 
volume limit. 
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Figure 5: The loop order attainable with a prescribed fraction a = -l(x), .2(+), .5 (*) compared 
with the estimate £ S( i o. 
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